Viral reprogramming of host transcription initiation.

Viruses are master remodelers of the host cell environment in support of infection and virus production. For example, viruses typically regulate cell gene expression through modulating canonical cell promoter activity. Here, we show that Epstein Barr virus (EBV) replication causes 'de novo' transcription initiation at 29674 new transcription start sites throughout the cell genome. De novo transcription initiation is facilitated in part by the unique properties of the viral pre-initiation complex (vPIC) that binds a TATT[T/A]AA, TATA box-like sequence and activates transcription with minimal support by additional transcription factors. Other de novo promoters are driven by the viral transcription factors, Zta and Rta and are influenced by directional proximity to existing canonical cell promoters, a configuration that fosters transcription through existing promoters and transcriptional interference. These studies reveal a new way that viruses interact with the host transcriptome to inhibit host gene expression and they shed light on primal features driving eukaryotic promoter function.

Alterations in gene expression and microbiome composition upon calcium-sensing receptor deletion in the mouse esophagus.

The calcium-sensing receptor (CaSR), a G protein-coupled receptor, regulates Ca2+ concentration in plasma by regulating parathyroid hormone secretion. In other tissues, it is reported to play roles in cellular differentiation and migration and in secretion and absorption. We reported previously that CaSR can be conditionally deleted in the mouse esophagus. This conditional knockout (KO) (EsoCaSR-/-) model showed a significant reduction in the levels of adherens and tight junction proteins and had a marked buildup of bacteria on the luminal esophageal surface. To further examine the role of CaSR, we used RNA sequencing to determine gene expression profiles in esophageal epithelia of control and EsoCaSR-/-mice RNA Seq data indicated upregulation of gene sets involved in DNA replication and cell cycle in EsoCaSR-/-. This is accompanied by the downregulation of gene sets involved in the innate immune response and protein homeostasis including peptide elongation and protein trafficking. Ingenuity pathway analysis (IPA) demonstrated that these genes are mapped to important biological networks including calcium and Ras homologus A (RhoA) signaling pathways. To further explore the bacterial buildup in EsoCaSR-/- esophageal tissue, 16S sequencing of the mucosal-associated bacterial microbiome was performed. Three bacterial species, g_Rodentibacter, s_Rodentibacter_unclassified, and s_Lactobacillus_hilgardi were significantly increased in EsoCaSR-/-. Furthermore, metagenomic analysis of 16S sequences indicated that pathways related to oxidative phosphorylation and metabolism were downregulated in EsoCaSR-/- tissues. These data demonstrate that CaSR impacts major pathways of cell proliferation, differentiation, cell cycle, and innate immune response in esophageal epithelium. The disruption of these pathways causes inflammation and significant modifications of the microbiome.NEW & NOTEWORTHY Calcium-sensing receptor (CaSR) plays a significant role in maintaining the barrier function of esophageal epithelium. Using RNA sequencing, we show that conditional deletion of CaSR from mouse esophagus causes upregulation of genes involved in DNA replication and cell cycle and downregulation of genes involved in the innate immune response, protein translation, and cellular protein synthesis. Pathway analysis shows disruption of signaling pathways of calcium and actin cytoskeleton. These changes caused inflammation and esophageal dysbiosis.

SpliceTools, a suite of downstream RNA splicing analysis tools to investigate mechanisms and impact of alternative splicing.

As a fundamental aspect of normal cell signaling and disease states, there is great interest in determining alternative splicing (AS) changes in physiologic, pathologic, and pharmacologic settings. High throughput RNA sequencing and specialized software to detect AS has greatly enhanced our ability to determine transcriptome-wide splicing changes. Despite the richness of this data, deriving meaning from sometimes thousands of AS events is a substantial bottleneck for most investigators. We present SpliceTools, a suite of data processing modules that arms investigators with the ability to quickly produce summary statistics, mechanistic insights, and functional significance of AS changes through command line or through an online user interface. Utilizing RNA-seq datasets for 186 RNA binding protein knockdowns, nonsense mediated RNA decay inhibition, and pharmacologic splicing inhibition, we illustrate the utility of SpliceTools to distinguish splicing disruption from regulated transcript isoform changes, we show the broad transcriptome footprint of the pharmacologic splicing inhibitor, indisulam, we illustrate the utility in uncovering mechanistic underpinnings of splicing inhibition, we identify predicted neo-epitopes in pharmacologic splicing inhibition, and we show the impact of splicing alterations induced by indisulam on cell cycle progression. Together, SpliceTools puts rapid and easy downstream analysis at the fingertips of any investigator studying AS.

Breast cancer cells survive chemotherapy by activating targetable immune-modulatory programs characterized by PD-L1 or CD80.

Breast cancer cells must avoid intrinsic and extrinsic cell death to relapse following chemotherapy. Entering senescence enables survival from mitotic catastrophe, apoptosis and nutrient deprivation, but mechanisms of immune evasion are poorly understood. Here we show that breast tumors surviving chemotherapy activate complex programs of immune modulation. Characterization of residual disease revealed distinct tumor cell populations. The first population was characterized by interferon response genes, typified by Cd274, whose expression required chemotherapy to enhance chromatin accessibility, enabling recruitment of IRF1 transcription factor. A second population was characterized by p53 signaling, typified by CD80 expression. Treating mammary tumors with chemotherapy followed by targeting the PD-L1 and/or CD80 axes resulted in marked accumulation of T cells and improved response; however, even combination strategies failed to fully eradicate tumors in the majority of cases. Our findings reveal the challenge of eliminating residual disease populated by senescent cells expressing redundant immune inhibitory pathways and highlight the need for rational immune targeting strategies.

Ebola virus delta peptide is an enterotoxin.

During the 2013-2016 West African (WA) Ebola virus (EBOV) outbreak, severe gastrointestinal symptoms were common in patients and associated with poor outcome. Delta peptide is a conserved product of post-translational processing of the abundant EBOV soluble glycoprotein (sGP). The murine ligated ileal loop model was used to demonstrate that delta peptide is a potent enterotoxin. Dramatic intestinal fluid accumulation follows injection of biologically relevant amounts of delta peptide into ileal loops, along with gross alteration of villous architecture and loss of goblet cells. Transcriptomic analyses show that delta peptide triggers damage response and cell survival pathways and downregulates expression of transporters and exchangers. Induction of diarrhea by delta peptide occurs via cellular damage and regulation of genes that encode proteins involved in fluid secretion. While distinct differences exist between the ileal loop murine model and EBOV infection in humans, these results suggest that delta peptide may contribute to EBOV-induced gastrointestinal pathology.

Reversal of splicing infidelity is a pre-activation step in B cell differentiation.

Introduction: B cell activation and differentiation is central to the adaptive immune response. Changes in exon usage can have major impacts on cellular signaling and differentiation but have not been systematically explored in differentiating B cells. Methods: We analyzed exon usage and intron retention in RNA-Seq data from subsets of human B cells at various stages of differentiation, and in an in vitro laboratory model of B cell activation and differentiation (Epstein Barr virus infection). Results: Blood naïve B cells were found to have an unusual splicing profile, with unannotated splicing events in over 30% of expressed genes. Splicing changed substantially upon naïve B cell entry into secondary lymphoid tissue and before activation, involving significant increases in exon commitment and reductions in intron retention. These changes preferentially involved short introns with weak splice sites and were likely mediated by an overall increase in splicing efficiency induced by the lymphoid environment. The majority of transcripts affected by splicing changes showed restoration of encoded conserved protein domains and/or reduced targeting to the nonsense-mediated decay pathway. Affected genes were enriched in functionally important immune cell activation pathways such as antigen-mediated signaling, cell cycle control and mRNA processing and splicing. Discussion: Functional observations from donor B cell subsets in progressive states of differentiation and from timecourse experiments using the in vitro model suggest that these widespread changes in mRNA splicing play a role in preparing naïve B cells for the decisive step of antigen-mediated activation and differentiation.

Cancers from Novel Pole-Mutant Mouse Models Provide Insights into Polymerase-Mediated Hypermutagenesis and Immune Checkpoint Blockade.

POLE mutations are a major cause of hypermutant cancers, yet questions remain regarding mechanisms of tumorigenesis, genotype-phenotype correlation, and therapeutic considerations. In this study, we establish mouse models harboring cancer-associated POLE mutations P286R and S459F, which cause rapid albeit distinct time to cancer initiation in vivo, independent of their exonuclease activity. Mouse and human correlates enabled novel stratification of POLE mutations into three groups based on clinical phenotype and mutagenicity. Cancers driven by these mutations displayed striking resemblance to the human ultrahypermutation and specific signatures. Furthermore, Pole-driven cancers exhibited a continuous and stochastic mutagenesis mechanism, resulting in intertumoral and intratumoral heterogeneity. Checkpoint blockade did not prevent Pole lymphomas, but rather likely promoted lymphomagenesis as observed in humans. These observations provide insights into the carcinogenesis of POLE-driven tumors and valuable information for genetic counseling, surveillance, and immunotherapy for patients. SIGNIFICANCE: Two mouse models of polymerase exonuclease deficiency shed light on mechanisms of mutation accumulation and considerations for immunotherapy.See related commentary by Wisdom and Kirsch p. 5459.

BH3 mimetics selectively eliminate chemotherapy-induced senescent cells and improve response in TP53 wild-type breast cancer.

TP53 wild-type breast tumors rarely undergo a complete pathological response after chemotherapy treatment. These patients have an extremely poor survival rate and studies show these tumors preferentially undergo senescence instead of apoptosis. These senescent cells persist after chemotherapy and secrete cytokines and chemokines comprising the senescence associated secretory phenotype, which promotes survival, proliferation, and metastasis. We hypothesized that eliminating senescent tumor cells would improve chemotherapy response and extend survival. Previous studies have shown "senolytic" agents selectively kill senescent normal cells, but their efficacy in killing chemotherapy-induced senescent cancer cells is unknown. We show that ABT-263, a BH3 mimetic that targets antiapoptotic proteins BCL2/BCL-XL/BCL-W, had no effect on proliferating cells, but rapidly and selectively induced apoptosis in a subset of chemotherapy-treated cancer cells, though sensitivity required days to develop. Low NOXA expression conferred resistance to ABT-263 in some cells, necessitating additional MCL1 inhibition. Gene editing confirmed breast cancer cells relied on BCL-XL or BCL-XL/MCL1 for survival in senescence. In a mouse model of breast cancer, ABT-263 treatment following chemotherapy led to apoptosis, greater tumor regression, and longer survival. Our results reveal cancer cells that have survived chemotherapy by entering senescence can be eliminated using BH3 mimetic drugs that target BCL-XL or BCL-XL/MCL1. These drugs could help minimize residual disease and extend survival in breast cancer patients that otherwise have a poor prognosis and are most in need of improved therapies.

Chemotherapy-induced senescent cancer cells engulf other cells to enhance their survival.

In chemotherapy-treated breast cancer, wild-type p53 preferentially induces senescence over apoptosis, resulting in a persisting cell population constituting residual disease that drives relapse and poor patient survival via the senescence-associated secretory phenotype. Understanding the properties of tumor cells that allow survival after chemotherapy treatment is paramount. Using time-lapse and confocal microscopy to observe interactions of cells in treated tumors, we show here that chemotherapy-induced senescent cells frequently engulf both neighboring senescent or nonsenescent tumor cells at a remarkable frequency. Engulfed cells are processed through the lysosome and broken down, and cells that have engulfed others obtain a survival advantage. Gene expression analysis showed a marked up-regulation of conserved macrophage-like program of engulfment in chemotherapy-induced senescent cell lines and tumors. Our data suggest compelling explanations for how senescent cells persist in dormancy, how they manage the metabolically expensive process of cytokine production that drives relapse in those tumors that respond the worst, and a function for their expanded lysosomal compartment.

Gammaherpesvirus RNAs Come Full Circle.

After an adaptive immune response is mounted, gammaherpesviruses achieve persistence through the utilization of viral noncoding RNAs to craft a suitable host cell environment in an immunologically transparent manner. While gammaherpesvirus long noncoding RNAs (lncRNAs) and microRNAs have been recognized for some time and have been actively investigated, a recent spate of reports have now identified repertoires of the circular RNA (circRNA) class of noncoding RNAs in both the lymphocryptovirus and rhadinovirus genera of gammaherpesviruses. Despite the recent nature of these findings, the detection of circRNAs across viruses and viral gene expression programs, the conservation of some viral circRNAs, and their detection in the clinical setting already raises the spectrum of functional importance in gammaherpesvirus biology and associated malignancies. Here, we provide an overview of currently known gammaherpesvirus circular RNAs and discuss reported physical and contextual properties that may be germane to future functional studies. With the Epstein-Barr virus (EBV) circRNAome being the most extensively studied to date, our discussions will be weighted toward EBV circRNAs while also addressing circRNAs discovered in the rhesus macaque lymphocryptovirus (rLCV), the Kaposi's sarcoma herpesvirus (KSHV), and the murid gammaherpesvirus 68 (MHV68). We hope that this will help set the stage for future investigations into the functions and relevance of this new class of viral noncoding RNAs in infection and disease.

Comparative Analysis of Gammaherpesvirus Circular RNA Repertoires: Conserved and Unique Viral Circular RNAs.

Recent studies have identified circular RNAs (circRNAs) expressed from the Epstein-Barr virus (EBV) and Kaposi's sarcoma herpesvirus (KSHV) human DNA tumor viruses. To gain initial insights into the potential relevance of EBV circRNAs in virus biology and disease, we assessed the circRNAome of the interspecies homologue rhesus macaque lymphocryptovirus (rLCV) in a naturally occurring lymphoma from a simian immunodeficiency virus (SIV)-infected rhesus macaque. This analysis revealed rLCV orthologues of the latency-associated EBV circular RNAs circRPMS1_E4_E3a and circEBNA_U. Also identified in two samples displaying unusually high lytic gene expression was a novel rLCV circRNA that contains both conserved and rLCV-specific RPMS1 exons and whose backsplice junctions flank an rLCV lytic origin of replication (OriLyt). Analysis of a lytic infection model for the murid herpesvirus 68 (MHV68) rhadinovirus identified a cluster of circRNAs near an MHV68 lytic origin of replication, with the most abundant of these, circM11_ORF69, spanning the OriLyt. Lastly, analysis of KSHV latency and reactivation models revealed the latency associated circRNA originating from the vIRF4 gene as the predominant viral circRNA. Together, the results of this study broaden our appreciation for circRNA repertoires in the Lymphocryptovirus and Rhadinovirus genera of gammaherpesviruses and provide evolutionary support for viral circRNA functions in latency and viral replication.IMPORTANCE Infection with oncogenic gammaherpesviruses leads to long-term viral persistence through a dynamic interplay between the virus and the host immune system. Critical for remodeling of the host cell environment after the immune responses are viral noncoding RNAs that modulate host signaling pathways without attracting adaptive immune recognition. Despite the importance of noncoding RNAs in persistent infection, the circRNA class of noncoding RNAs has only recently been identified in gammaherpesviruses. Accordingly, their roles in virus infection and associated oncogenesis are unknown. Here we report evolutionary conservation of EBV-encoded circRNAs determined by assessing the circRNAome in rLCV-infected lymphomas from an SIV-infected rhesus macaque, and we report latent and lytic circRNAs from KSHV and MHV68. These experiments demonstrate utilization of the circular RNA class of RNAs across 4 members of the gammaherpesvirus subfamily, and they identify orthologues and potential homoplastic circRNAs, implying conserved circRNA functions in virus biology and associated malignancies.

Breast cancer survival predicted by TP53 mutation status differs markedly depending on treatment.

BACKGROUND: Previous studies on the role of TP53 mutation in breast cancer treatment response and survival are contradictory and inconclusive, limited by the use of different endpoints to determine clinical significance and by small sample sizes that prohibit stratification by treatment. METHODS: We utilized large datasets to examine overall survival according to TP53 mutation status in patients across multiple clinical features and treatments. RESULTS: Confirming other studies, we found that in all patients and in hormone therapy-treated patients, TP53 wild-type status conferred superior 5-year overall survival, but survival curves crossed at 10 or more years. In contrast, further stratification within the large dataset revealed that in patients receiving chemotherapy and no hormone therapy, wild-type TP53 status conferred remarkably poor overall survival. This previously unrecognized inferior survival is consistent with p53 inducing arrest/senescence instead of apoptosis. Addition of hormone therapy to chemotherapy improved survival notably in patients with TP53 wild-type tumors, but not mutant, suggesting hormone therapy could eradicate arrested/senescent cells. Testing this, we found that estrogen receptor-positive, TP53 wild-type breast cancer cells that were made senescent by doxorubicin treatment were sensitive to tamoxifen. CONCLUSIONS: The poor survival of chemotherapy-treated patients with TP53 wild-type tumors may be improved by strategies to eliminate senescent cells, including the addition of hormone therapy when appropriate.

p53 Mediates Vast Gene Expression Changes That Contribute to Poor Chemotherapeutic Response in a Mouse Model of Breast Cancer.

p53 is a transcription factor that regulates expression of genes involved in cell cycle arrest, senescence, and apoptosis. TP53 harbors mutations that inactivate its transcriptional activity in roughly 30% of breast cancers, and these tumors are much more likely to undergo a pathological complete response to chemotherapy. Thus, the gene expression program activated by wild-type p53 contributes to a poor response. We used an in vivo genetic model system to comprehensively define the p53- and p21-dependent genes and pathways modulated in tumors following doxorubicin treatment. We identified genes differentially expressed in spontaneous mammary tumors harvested from treated MMTV-Wnt1 mice that respond poorly (Trp53+/+) or favorably (Trp53-null) and those that lack the critical senescence/arrest p53 target gene Cdkn1a. Trp53 wild-type tumors differentially expressed nearly 10-fold more genes than Trp53-null tumors after treatment. Pathway analyses showed that genes involved in cell cycle, senescence, and inflammation were enriched in treated Trp53 wild-type tumors; however, no genes/pathways were identified that adequately explain the superior cell death/tumor regression observed in Trp53-null tumors. Cdkn1a-null tumors that retained arrest capacity (responded poorly) and those that proliferated (responded well) after treatment had remarkably different gene regulation. For instance, Cdkn1a-null tumors that arrested upregulated Cdkn2a (p16), suggesting an alternative, p21-independent route to arrest. Live animal imaging of longitudinal gene expression of a senescence/inflammation gene reporter in Trp53+/+ tumors showed induction during and after chemotherapy treatment, while tumors were arrested, but expression rapidly diminished immediately upon relapse.

Increased activin bioavailability enhances hepatic insulin sensitivity while inducing hepatic steatosis in male mice.

The development of insulin resistance is tightly linked to fatty liver disease and is considered a major health concern worldwide, although their mechanistic relationship remains controversial. Activin has emerging roles in nutrient homeostasis, but its metabolic effects on hepatocytes remain unknown. In this study, we investigated the effects of increased endogenous activin bioactivity on hepatic nutrient homeostasis by creating mice with inactivating mutations that deplete the circulating activin antagonists follistatin-like-3 (FSTL3) or the follistatin 315 isoform (FST315; FST288-only mice). We investigated liver histology and lipid content, hepatic insulin sensitivity, and metabolic gene expression including the HepG2 cell and primary hepatocyte response to activin treatment. Both FSTL3-knockout and FST288-only mice had extensive hepatic steatosis and elevated hepatic triglyceride content. Unexpectedly, insulin signaling, as assessed by phospho-Akt (a.k.a. protein kinase B), was enhanced in both mouse models. Pretreatment of HepG2 cells with activin A increased their response to subsequent insulin challenge. Gene expression analysis suggests that increased lipid uptake, enhanced de novo lipid synthesis, decreased lipolysis, and/or enhanced glucose uptake contribute to increased hepatic triglyceride content in these models. However, activin treatment recapitulated only some of these gene changes, suggesting that increased activin bioactivity may be only partially responsible for this phenotype. Nevertheless, our results indicate that activin enhances hepatocyte insulin response, which ultimately leads to hepatic steatosis despite the increased insulin sensitivity. Thus, regulation of activin bioactivity is critical for maintaining normal liver lipid homeostasis and response to insulin, whereas activin agonists may be useful for increasing liver insulin sensitivity.

Differential synthesis and action of TGFß superfamily ligands in mouse and rat islets.

Members of the TGFß superfamily, including activins and TGFß, modulate glucose-stimulated insulin secretion (GSIS) in vitro using rat islets while genetic manipulations that reduce TGFß superfamily signaling in vivo in mice produced hypoplastic islets and/or hyperglycemia. Moreover, deletion of Fstl3, an antagonist of activin and myostatin, resulted in enlarged islets and ß-cell hyperplasia. These studies suggest that endogenous TGFß superfamily ligands regulate ß-cell generation and/or function. To test this hypothesis, we examined endogenous TGFß ligand synthesis and action in isolated rat and mouse islets. We found that activin A, TGFß1, and myostatin treatment enhanced rat islet GSIS but none of the ligands tested enhanced GSIS in mouse islets. However, follistatin inhibited GSIS, consistent with a role for endogenous TGFß superfamily ligands in regulating insulin secretion. Endogenous expression of TGFß superfamily members was different in rat and mouse islets with myostatin being highly expressed in mouse islets and not detectable in rats. These results indicate that TGFß superfamily members directly regulate islet function in a species-specific manner while the ligands produced by islets differ between mice and rats. The lack of in vitro actions of ligands on mouse islets may be mechanical or result from species-specific actions of these ligands.